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Abstract 

Locally, isoperimetric problems on Riemannian surfaces are sub-Riemannian 
problems in dimension 3. The particular case of Dido problems corresponds to 
a class of singular contact sub-Riemannian metrics : metrics which have the chara- 
teristic vertor field as symmetry. We give a classification of the generic conjugate 
loci (i.e. classification of generic singularities of the exponential mapping) of a 
1-parameter family of 3-d contact sub-Riemannian metrics associated to a 1-parameter 
family of Dido Riemannian problems. 



1 Introduction 

1.1 Sub-Riemannian metrics under consideration 

Let M. be any smooth 3- dimensional manifold. Let TA4 and T*M. be respectively its 
tangent and cotangent bundle and let 7r : T*j\4 — > M. be the natural projection. 
A contact sub-Riemannian metric on M. is a couple S = (A,g), where 

A = {A q } qeM ,A q C T q M, 



is a contact structure on M. and 

9 ■ A -> R +1 

a Riemannian metric on A. 

Since A is a contact structure and therefore nonintegrable, (A, g) defines a distance d on 

M !. 

A contact sub-Riemannian metric being given on Ai, there is a 1-form uj and a vector 
field v on A4, both determined up to orientation by the conditions: 



1. The distribution A is the kernel of uj (A = ker uj), 

2. the 2-form du restricted to A is the volume form V, 
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3. the vector field Z is such that u(Z) = 1, i(Z) du = 0. 

The vector field Z is the characteristic vector field of the contact sub-Riemannian struc- 
ture. 

The cotangent bundle has the standard symplectic structure. Consider the Hamiltonian 
TC of the metric. It is defined as follows: associated to (A, g), there is a cometric on 
T*M : 

H(il>) = I sup ( ) , where i(,(q) G T*M. 

2„gA\{o} \ \\u\\ g J 

On fibers of T*M, H is a positive semi-definite quadratic form, the kernel of which is the 
annihilator of A. 

In the remainder of this paper we will assume that: 

1. The manifold M. is an open subset of M 3 , containing the origin 0, 

2. the distribution A is specified by an orthonormal frame field (F, G) of the metric 
g, where F and G are two vector fields defined on Ai, 

3. for any q G A4, F(q), G(q), and [F,G](q) are linearly independent (the contact 
condition). 

Hence the Hamiltonian H is given by: 

n= l -{^{F)f + ^{G)f). 

In this paper, we deal with contact sub-Riemannian structures that are invariant under 
the action of the one parameter group generated by the characteristic vector field Z ( i.e. 
sub-Riemannian structures having a symmetry Z). The reasons for doing so will become 
clear when we show that: to characterize minimizing geodesies for this class of metrics 
is a dual sub-Riemannian reformulation of the classical isoperimetric problem of Dido on 
Riemannian surfaces. 

Let (M, g) be a Riemannian surface. We consider on M the following iso-area problem: 
First we fix two points z Ql z\ G M and a curve 

I: [0,l]->M,r(0) = ;zo,Z:(l) = Zi. 

We are then faced with the following question: 
Which curves 

t :[0,l]-M, 
connecting z and z 1 , such that the area 

A= ( V 
Jn 

of the domain Q encircled by I and i is prescribed, have minimal Riemannian length? 
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Denoting such an iso-area problem (M, g, V), it is stated (see for instance [[|) that 
this problem can be reformulated locally in terms of three dimensional sub-Riemannian 
geometry. 

More precisely, we can associate a germ of Dido Problems {M,g, V) ZQ with a germ of 
an oriented sub-Riemannian structure on M. with a symmetry Z, denoted by: 

(M,A,g,Z) qo , 

where, 

M = M x R and q = {z , 0} 

d 

Let V = da and Z = — — , then 

ow 

A = ker(dw + a). 

If 

S:[Q,1]^M 
is an admissible curve for A with fixed endpoints 

5(0) = q and 5(1) = q x = (^i,Wi), 

hence: 

The sub-Riemannian length of 5 is the Riemannian length of its projection 5* on M and 

1 ,d8* 



Wi = a(——)dT. 
Jo dr 



Therefore projections of sub-Riemannian length minimizers are solutions of the iso-area 
problem (M, g, V). 

The sub-Riemannian structures associated with (M, g, V)) or simply (M, g) are called the 
isos-R-structure and are denoted by IsosR. 

The main aim of this paper is to classify generic conjugate loci for IsosR-metrics. 

For some reasons which will become clear in the next section we will in fact consider 
1-parameter families of Dido Riemannian problems: 

A G / C R i-> (M,g(X)), I an interval. 

We will denote the families of Jsosi?-metrics associated with (M,g(X)) by F-IsosR. 

1.2 Normal forms for Riemannian metrics on surfaces 

Recall the following result. 
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Theorem 1.1 Let (M,g) be a Riemannian surface. The metric g has the following nor- 
mal form. In normal coordinates with pole 0, there is an orthonormal frame (unique up 
to rotations and up to the action of SO (2) in T M ), (X, Y): 

(SNF) I 

Where (3(x,y) is a smooth function . Moreover 

/3(0,0) = &«, = ifco. 

o 

The constant ko being the gaussian curvature of M at 

Remark 1.1 In the normal coordinates specified above, the h th - differentials (3 k = D k /3 | o 
are homogeneous polynomials of degree k in x,y and define symmetric covariant tensors 
of degree k on TqM. 

These tensors, denoted /3&, do not depend on the orientation and are invariants of the 
Riemannian structure. 

Decomposition of tensor fields Let S k T*M denote the bundle of covariant sym- 
metric tensors of degree k on TM. Due to the action of SO (2) on the typical fiber of TM, 
we have the following decomposition of S k T*M into isotypic components @|: 

S k T*M = © A/ (5 fc T*M) i; 
where (S k T*M)j is the component relative to the j th power of the basic character e iV , 

(< = >/=!). 

Therefore, if (3 n e S k T*M then: 

P n = ^Pru, where G (S k T*M) r 
i 

1.3 Lemmas 

Let M be any open subset of M 2 . Let O(M) and J n O(M) be respectively the bundle 
of orthonormal frames of Riemannian metrics on M and the vector bundle of n-jets of 
elements of O(M). 

We denote by IT, : 0(M) M and IT : J n O{M) -> M the standard projections. 
Let 

B n = ®{S k T*M | < k < n}, B' 1 = {0} 
For n > 0, let us define the following map: 

rr : j n+2 o(M) -> B n -\ 



X dy)- 
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by its restriction to the fiber IT" 1 ^) of the bundle IT : J n+2 0(M) -> M: 

where /3 n_1 is the representative of the tensors /? n _i 90 in the unique normal coordinates 
system at q such that: 

(X(q ),Y(q )) = (X(q ),Y(qo)), 
where (X, Y) is the normal form of (X, Y) at go- 
Lemma 1.1 If we fix (X,Y) E 0(M), then: 

U n (X,Y) : M -> B n -\ U n (X,Y)(q ) = U n (j^ 2 (X),j^ 2 (Y)). 
is a surjective submersion. 



The proof of lemma |LT] is the result of computations not given here. 
Definition 1.1 A 1-parameter family T of orthonormal frames is given by: 

{{X x ,Y x )eO(M),\eI}. 

Lemma 1.2 If we fix (X, Y) E J 7 , 

fl n (X,Y) : M x / - B n -\ U n (X,Y)(q ) = U n (j^ o) (X) , j^jY)) 

is a surjective submersion. 



Lemma |1.2| is a consequence of lemma [0| , and of the fact that J n T submerses into 
J n O(M). 

The algebraic lemmas (lemma IT and lemma 1.2) are the main elements in establishing 
the genericity results here. 

Since we have two lemmas, we have the following two choices: We can envisage an isolated 
Dido problem and use lemma |1 . 1| , or a family of such problems depending on 1-parameter 



by using lemma 1.2 



It is obvious that there are nongeneric situations in the first case, which become generic 
in the second. Naturally we will consider 1-parameter families. 

1.4 Exponential mapping 

Let us first recall briefly the previous results of M, M and |T||, which concern generic 
contact sub-Riemannian structures in dimension 3. 
The following result is stated in M. 
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Theorem 1.2 Consider a germ at of contact sub-Riemannian metric (A,g). There is 
up to orientation, and up to the action of SO (2) on A(0), a unique coordinate system 
(normal coordinates) with respect to which the metric has an orthonormal frame (F,G) 
of the following form: 

\ ^ = A + ^(^A-^) + K 1 + 7)i) 

(NF) \ (2) 
[ G = 4-i(^(y4-xJ) + l(l+7)A). 

Where (3 and 7 are smooth functions ; satisfying the following boundary conditions: 

,3(0, 0, w) = 7(0, 0, w) = ||(0, 0, u>) = |?(0, 0, w) = 0. 

Indeed, {NF) can be obtained even when A is not a contact structure. In fact the normal 
coordinates are coordinates w.r.t. local coordinates charts 

(U r , (w,p,q) o (yroexp H)' 1 ), 

where (p, q, r) are dual coordinates in T*A4 and T is a smooth curve transversal to A. In 
the contact case, such a curve T can be taken as the integral curve of Z. 

In T-normal coordinates, geodesies through T, satisfying transversality conditions 
w.r.t. T are straight lines contained in the plane {w = c}, where c is a constant ([![]). The 
set of points x of A4, that are at a distance s of V, form a smooth cylinder for small 
enough values of s. 



1.4.1 Approximation of generic conjugate loci and stability 

In p| and |J, it is proved that for contact sub-Riemannian metrics, there are in essence 
two generic situations. In the two cases we have a representation of the exponential 
mapping as a suspension of a mapping between surfaces. 

The first case happens on a complement of a smooth possibly empty curve C of M.. In 
this instance, the 3-jet (exp 3 ) is a sufficient jet of the exponential mapping. 
This result follows the fact that the suspension of (exp 3 ) is a "Whitney map": that is a 
stable mapping in the sense of Thorn-Mather ( |T3 ) . 

Hence the 3-jet of the conjugate loci is sufficient to describe it and the exponential mapping 
is determined by the 3-jet of the metric, in a neighborhood of its singular locus. 
The situation on the curve C is more complex. The full conjugate loci CL = CL, splits 
into two semi-conjugate loci CL + and CL~ corresponding to w > or w < in the 
normal coordinates. 

The intersection CL^ of CL^ with the planes {w = c 7^ 0} for small enough values of 
\w\, is a closed curve with 6 cusps and self-intersections. Hence the suspension of the 
exponential mapping is not a "Whitney map" . 

To conclude that the exponential mapping approximation is stable, it is necessary to show 
that all generic self-intersections are transversal. 

In order to classify the semi-conjugate loci, the authors of define the symbol S of CL^ 
in the following manner. 

A symbol S is a sequence of 6 rational numbers, S — (s%, . . . , Sq) modulo cyclic permuta- 
tion and reflection. If we follow the curve 5 starting from a cusp point: Si gives half the 
number of self-intersections between the i th and the (i + l) th cusp point. 
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Theorem 1.3 1) At generic points of the curve C, the possible symbols for generic semi- 
conjugate loci are: 



Si = (2, 1,1, 2, 1,0), S 2 = (2, 1,1, 1,1,1), 

S 3 = (0,1, 1,1, 1,1). 

2) There are two types of isolated points at C. 

a) At the first type the possible symbols are: 

si = S l = (2, 1, 1, 2, 1, 0), s; = S 2 = (2, 1,1,1, 1, 1), 

S* = S 3 = (0,1, 1,1, 1,1). 

But in that case , the exponential mapping is determined by the higher order jet of the 
metric (7-jet) than at generic points of C (5-jet). 

b) At the second type the possible symbols are: 

^4 = (|, |, 1, 0, 0, 1) S 5 = (1, 2, 2> 1? lj 1) 

S 6 = (§,§, 1,1, 0,1) ^=(2,1, |,2,0,0). 

However, this result, stated also in Q, only gives classification of possible semi-conjugate 
loci. The problem of classifying possible full conjugate loci is therefore not entirely solved. 



1.5 Normal Dido coordinates 

Theorem 1.4 Given (A,g) a germ at the origin of an isoperimetric sub-Riemannian 
metric. There is an unique coordinates system (normal Dido coordinate), up to orientation 
and up to the action of SO (2) on A(0), with respect to which the metric has an orthonormal 
frame (Fj, Gi) of the following form: 

{ ^ = X +§7£ 
(INF) { (3) 

Where : (X, Y) is the normal form of the Riemannian metric on M and 

? = < 1 + <^W lWT7b<^ - (4) 



This theorem is not proved here. Observe that in comparison with theorem |1.2| , the only 
work we have to do for this, is to slightly modify the T-normal coordinates. This done, 
checking formula ^ is very easy. 

In accordance with |5| the jF-Jsos-R-structures are nongeneric. However this is fairly 
interesting that as it is in the case of generic contact sub-Riemannian metrics, our work 
follows almost the same steps and we prove similar results. 
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1.6 Statement of our main results and outline of the paper 

Our work is organized as follows. In the second section, we give the notations for the 
main invariants. We end this section by summarizing the basic properties of the relevant 



jet of the conjugate locus computed in appendix p73 . 

In section 3, we state our genericity results. We highlight two principal invariants denoted 
by r 2 and r 3 . More precisely we prove : 

For a generic element £ of T (for the Whitney topology), the set of points of M x I on 
which the invariant r 2 vanishes is a smooth possibly empty curve C . 
On the curve C the invariant r 3 does not vanish. 



Using higher order invariants, theorems |3.2| and |3.3| show that the curve C is partitioned 
into two subsets: a discrete subset (isolated points of C) and its complement (generic 
point of C). 

In section 4, we state our stability results (see theorems |0, [4.2| , |4.3| and |4.4j ). We can 
summarize these results as follow: 

If r 2 is nonzero, the exponential mapping is determined by the 5-jet of the metric, in a 



neighborhood of its singular locus. . This is the purpose of theorem p~T . 

On C, the situation is by far more complex and particularly delicate at isolated points on 

this curve. 



For the sake of clarity theorem [4.2| deals only with generic points on the curve C. Thus 
theorems O and O are devoted to isolated points on C. 

It is clear that for J-'-IsosR, we can classify the semi conjugate loci as in theorem |1.3| , but 
we need higher order jets. 

What is more, for F-IsosR, our symbols are the same for CL + and CL~ . It follows that 
we can classify the full conjugate loci. 

In fact we state the following results. 

Theorem 1.5 1) At generic points of the curve C , the possible symbols for generic semi- 
conjugate loci are: 

Si = (2, 1,1, 2, 1,0), £, = (2,1,1,1,1,1), 

S 3 = (0,1, 1,1, 1,1). 

2) There are two types of isolated points at C. 

a) At the first type the possible symbols are: 

S* = Si = (2, 1, 1, 2, 1, 0), S 2 * = S 2 = (2, 1, 1, 1, 1, 1), 

S 3 * = S 3 = (0,1, 1,1, 1,1). 

But in that case , the exponential mapping is determined by higher order jets of the 
metric (9-jet) than at generic points of C (7 -jet). 

b) At the second type the possible symbols are: 

S 4 = (|, |, 1, 0, 0, 1) S 5 = (1, |, |, 1, 1, 1) 
S 6 = (§, |,1, 1,0,1) S 7 = (2, |, |,2,0,0). 
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Theorem 1.6 For a fixed element of F-IsosR the possible symbols for generic conjugate 
loci are: 



Si = (2, 1,1, 2, 1,0), S 2 = (2, 1,1, 1,1,1), 

S 3 = (0,1, 1,1, 1,1). 

The section 5 is devoted to technical computations of the exponential mapping and 
the conjugate loci. 

2 Conjugate loci 

Let E = (A,g)o be a germ at the origin of an isoperimetric sub-Riemannian metric. We 
can restrict M. in order to obtain a neighborhood N of xq = (0,0,0), with respect to 
which £ is under normal form (INF). We denote by (p, q, r) the dual coordinate in T*N. 

The set C| = ir^i^o) H H~ l (\) is therefor the cylinder {p 2 + q 2 = 1}. 

Let expsH : T*IR 3 — > T*M 3 , denotes the Hamiltonian flow of H at time s. If N is small 

enough, point of the conjugate loci in iV appear only along geodesies {n o exp sH(X X0 )} 

i 

for X xo in a certain neighborhood of infinity TV* C C£ 0. 



In appendix ^]3| we compute relevant jets of CL in a suitable coordinates system (h, ip), 



with respect to which: 

CL : iV* -> N 
has the following expansion with respect to h: 

9 

CL( V , h) = (x(if, h), y(<p, h),h) = (J2 fifr) V + 0(h w ), h), 

i=4 

where /i = v/^, e = ±1 according to the sign of w, and (x, y, w) the normal coordinates. 

CL = CL + U CL- 



As we will see, for small enough of \w\, the intersection of CL* 1 with planes {w = c ^ 0} 
will be a close curve. In the more degenerate generic cases this curve will have self- 
intersections denoted by Isoself. 

2.1 Notation 

The following tensors have the leading part in our study. 

Let Pi = /?2 = (32,2 + ^2,0, and (3% = /%,3+/?3,i. Thus we have in our normal coordinates 
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Where : 



= R e {n{dx + idy)). 

02.2 = R e (r 2 (dx + idy) 2 ), (3 2 ,o = t (dx 2 + dy 2 ). 

03.3 = Re(r 3 (dx + idy) 3 ), (3 3A = R e (v (dx + idy))(dx 2 + dy 2 ). 

r\ = | r*i j (cos 0i, sin 9\ ) = |ri| (cos^ + zsin^i) 
r 2 = \r 2 \ (cos6> 2 , sin^) = \r 2 \ (cos6 2 + isin9 2 ) 
r 3 = \r 3 \ (sin ^3, -cos6» 3 ) = |r 3 | (sin6> 3 - icos# 3 ) 
v = —v\ +iv 2 . 



The remaining tensors are given in appendix 5.3 



2.2 Basic properties of f{ (B.P.) 

From now on, we summarize basic properties of the map fa when e = 1. For details see 
appendix FT3] 



Pi) r 2 + 0. 



1. fify) = 37r |ri| (— sin 9, cos0), 

2. / 4 is independant of <p, 

3. f 5 ((p) = 57r \r 2 \ (3 cos(v2 - 9 2 ) + cos(3<y? - 9 2 ), 3 sin(^ - 9 2 ) - sin(3<y3 - 9 2 )). 



P 2 ) r 2 = 0. 



1. / 4 (^) = 37r|r 1 |(-sin^ 1 ,cos^ 1 ), / B (p) = (0,0), 
2. 

/e<» = f (-25V2 + 90 |r 3 | cos(2</? + 9 3 ) + 45 |r 3 | cos(4<^ + 9 3 ) + 31b |n| sin 0i, 
-25^i - 90 |r 3 | sin(2<^ + 3 ) + 45 |r 3 | sin(4^ + 6> 3 ) - 316 N cos0i), 

3. f 6 ((p + 7r) = f 6 (<p), 

4. f 7 (ip + Ti) = -f 7 (<p) t 
5. 

^ A ^ = 12960tt 3 |r 3 | sin(3^ + 6> 3 ) |r?| sin2(^ - B x ), 
dtp dip 
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6. 

<!J 6 A f 7 = 1080tt 2 |r 3 | sin(3y? + 9 3 )P{cp) = V(<p), 



dip 

7. f 8 (tp + n) = f 8 ((f), 

8. f 9 {v + 7i) = -f 9 (<p), 
9. 



10. 



^A^ = \r 3 \ S m(3<p + 9 3 )(Q Q (<p)), 
dip dip 



-j- A /g = |r 3 | sin(3^ + G 3 ) (Q l (i P ) + Q 2 (<p))). 
dip 



Remark 2.1 P(<p) = A cos(2ip)+B sin(2y3)+C cos(4^) + D sin(4y?), where A, B, C, D 
are linear combinations of the coefficient of the invariants (3^ and r\. 
Contrary to the generic contact case, P(<p) is indepent on e (see appendix \5.ty ). 



Remark 2.2 The Qi's do not depend on the same invariants. 

In particular only Q 2 is dpendent on the coefficients of /3g. 
Let : 

A-iB C-iD 
2 ' P 2 

Let 

P((p) = ue 2 ^ + Ve~ 2t ^ + /ie 41 ^ + 7Ze" 4 ^. 
and 

P(z) = fiz 4 + h>z 3 + T>z + JJ 
Tc(tp) = |r 3 | sin(3v9 + 6> 3 ) 
and 

T(z) = r 3 z 3 + f 3 . 

In the remainder of this paper we will be concerned with the roots on the unit circle of 
these trigonometric polynomials. 
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3 Genericity results 



In the next section we will prove the following results : 

Theorem 3.1 For a generic element of J 7 (for the Whitney topology), the set on which 
02,2 vanishes is a smooth curve C ofM, 3 . 
Along C, /3 3)3 and fi do not vanish. 

Theorem 3.2 On the complement of a discrete subset of C : 

i ) P(z) has either two or four simple roots on the unit circle, 

ii ) P(z) and T(z) have no common roots on the unit circle. 

Theorem 3.3 At isolated point ofC, either : 

i) P(z) has a double root on the unit circle, which is not triple, and which is not a 
root ofT(z). 

The other roots of P(z) are simple ; they are on the unit circle and are not roots of 
T(z). 
or, 

ii) P(z) and T(z) have one and only one common simple root on the unit circle. 
The other roots of P(z) are simple. 

3.1 Proof of our genericity theorems 
3.1.1 The bad set 

Definition 3.1 B is the subset of B 6 such that : 

1) /?2,2 = &,3 = 

B\ is the subset of B & such that : 

1) &,2 = /i = 

Definition 3.2 B 2 is the subset of B 6 such that : 

1) fh,2 = 

2) P(z) has a triple root on the unit circle. 
B 3 is the subset of B 6 such that : 

1) fo,2 =0, 

2) A root ofT(z) is a double root of P(z) on the unit circle. 
B± is the subset of B 6 such that : 

1) P2,2=0, 

2) T(z) and P(z) have two common roots on the unit circle. 
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Definition 3.3 B 5 is the subset of B 6 such that 

P{z) has a double root on the unit circle 
B 5 is the subset of B 6 such that 

1) P2,2 = 0, 

P(z) and T(z) have one and only one common simple root on the unit circle 
(i.e. Res(r3, //, v) =Resultant of P(z) and T(z) has a simple root). 

Remark 3.1 Since the integral ofP(z) over its periodis zero, then, if /i ^ 0, the following 
facts are easy to check. 

1) P(z) has a root on the unit circle. 

2) P(z) does not have two double roots on the unit circle. 

3) If P(z) has a double root on the unit circle, the two other roots of P(z) are on the 
unit circle. 

3.1.2 Estimate of the codimension of B^'s 

Some appropriate computations on T(z) and P(z) allow us to rewrite the B k s as follows: 



B 


= {r 2 


= 0} n {r 3 = 0} 




B 1 


= {r 2 


= o}n{/i = o} 




B 2 


= {r 2 


= 0} n {4/i 2 z/ + z/ 3 = 




B 3 


= {r 2 


= 0} n {27u 3 r 3 fl + ( 


[9r 3 - 4:fif 3 ) 3 = 0} 


B 4 


= {r 2 


= 0} n {z/r 3 - jj,r 3 = 




B 5 


= {r 2 


= 0} n {27 Re 2 (/iP 2 ) 


- (4|^| 2 - \U\ 2 ) 3 : 


B 6 


= {r 2 


= 0} n {Res(/x,^ r 3 ) 


= 0} 



(5) 



0} 



3.1.3 Proof of the genericity 

We consider the B^s defined from the B^s by 



B k = (n^)- 1 ^) 

Now let : 

B k = {(X,Y)\J 9 (X,Y)(hB k } 

If B = C\ k B k ; then B is an open dense set for the Whitney topology in T . 
Standard arguments from transversality theory (see for instance Goresky-MacPherson 
fTT|), allow us to conclude that : 
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V(X,Y)eB, (fL^(X,Y)- l )(B k ) 

are smooth submanifolds or Whitney-stratified set of M x I of the same codimension 
as the -Bfc's. 

Since Bk's. are at least codimension 4 (resp codimension 3) for < k < 4 (resp 5 < k < 6), 



this ends the proofs of theorems kO, kO.and 3.3. 



Remark 3.2 If we consider a fixed element T\ of J 7 , then: for a generic element of 
J-'xq (f or the Whitney topology), 02,2 vanishes at isolated points ofM. 2 . 



Remark 3.3 at the isolated points of C specified in the part (i) of the theorem \3.3( , we 
have : 

Q^O, ? = 0,1,2 
Qo + Q2 ^ 

Denoting £? 5 the subset of B 6 x S" 1 (where S" 1 is the unit circle) defined by the following 
equations : 

dP 

(32,2 = 0, P(cp) = 0, —{<p) = 0, Qi{tp) = (resp Q + Q 2 = 0). 
dip 



Taking into account the remark |2.2| , we show as in [y] that the projection B 5 of B 5 on 
B 6 is a semi-algebraic subset of codimension four. 

4 Stability results 

Let CL be the conjugate locus mapping for a generic element of J r -IsosR (for the Whitney 
topology) and S a neighborhood of the singular locus at the source. 

Theorem 4.1 On the complement of the smooth curve C (/?2,2 7^ 0). 

There is a neighborhood U of S H {0 < h < c} (for small enough values of c), such 
that : 

The 5-jet of CL : CL? = f^h 4 + f 5 h 5 is a sufficient jet of CL on U. 
The restriction of the exponential mapping exp \U is 5-determined in h. 
exp \ U is determined by the 5-jet of the metric. 



The main argument for the proof of the theorem |4.1| is the fact that the suspension of the 
exponential mapping is a "Whitney map" 

In fact the intersection CL^ between CL and the planes w = c (for small enough values of 
|w|) is a closed curve which has fold-points, cusp-points (four) and without self- intersection 
(see figures 1 and 2). 
Hence the mapping 

M 4 + fch 5 + 0(h 6 ) 

is R.L-equivalent to 

M 4 + M 5 
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Figures 




CL b when r 2 ^ 0. 




CLt (when r 2 ^ 0) 



Theorem 4.2 There is an open dense subset O, (complement of a discrete subset) ofC, 
such that : 

Isoself = Isoself c + Isoselfi 

— Isoself c is the union of three curves. 
—Isoselfi is the union of two or four curves. 

All the self-intersections are transversal and are not dependent on e. 
There is a neighborhood U of § fl {0 < h < c} (for small enough values of c), such 
that : 

The 7-jet of CL : CL 7 = f^h 4 + f 6 h 6 + f 7 h 7 is a sufficient jet of CL on U. 
The restriction of the exponential mapping exp \ U is 7 '-determined in h. 
exp \ U is determined by the 9- jet of the metric. 

At the isolated points of C we have : 

Theorem 4.3 If P{ip) has a double root on the circle unity or P(ip) = and fi{tp) = 0, 
then : 

—Isoself c is the union of three curves. 

—Isoselfi is the union of two or four curves. 

All the self-intersections are transversal and are dependent on e. 

There is a neighborhood U of § fl {0 < h < c} (for small enough values of c), such 
that : 

The 9-jet of CL : CL 9 = f 4 h 4 + f 6 h 6 + f 7 h 7 + f s h 8 + f 9 h 9 is a sufficient jet of CL on 

U. 

The restriction of the exponential mapping exp \ U is 9-determined in h. 
exp \ U is determined by the 9-jet of the metric 



Theorem 4.4 When there is collision between Isoself c and Isoselfi, then : 
—Isoself c is the union of two curves. 
—Isoselfi is the union of one or two curves. 

All the self-intersections are transversal and are not dependent on e. 
There is a neighborhood U of § fl {0 < h < c} (for small enough values of c), such 
that : 

The 7-jet of CL : CL 7 = f A h 4 + f 6 h 6 + f 7 h 7 is a sufficient jet of CL on U. 
The restriction of the exponential mapping exp \ U is 7 -determined in h. 
exp \ U is determined by the 7-jet of the metric 

4.1 Proof of the stability results along the curve C 

4.1.1 Definition and characterization of the self-intersections of CL 

Let Isoself denote the set of self-intersections of CL on C, roughly speaking the set of 
(h, ipi, ip2) such that : 

1) <Pl ^ ¥2, 

2) CL(h,ipi) = CL(h,<p 2 ). 
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We are interested in the germ of Isoself along {h = 0}. A value of f is said to be 
adherent to Isoself if (0, tp, ip ) lies in the closure of Isoself for some (p . The set of some 
ip is denoted by A-Isoself. 

We know that if tp G A- Isoself then ip = ip + 7t G A- Isoself and there is not an other 
possibility (0j). 

Lemma 4.1 A-Isoself C {95 1 ^(95) = 0}. 

Assuming that (h, <p, <p + 7r + 5) G A-Isoself, with h > 0, and 5 small. Then : 

7 7 

+ 7T + = £/,fo>)tf + ^(/i 8 ). (6) 

1=4 i=4 

Since /4 is independent of 9? and /s^) = on C, dividing |6] by h G we obtain : 

fe(<P + - fM + h (f 7 (<P + S + 7T)- f 7 (<p + rr)) + h {f 7 { V + rr) - f 7 (<p)) + 0(h 2 ) = 
Using (B.P.) we obtain : 

Safety) - 2hf 7 {tp) + 0(/i 2 ) + 0(hS) = 
at/? 

Therefore if (p G A-Isoself , then ^(<£>) A ^(ty?) = ; and lemma [O] is proved. 

Remark 4.1 ^4s a consequence of this lemma A-Isoself splits into two subsets : 

A-Isoself = A-Isoself c U A-Isoselfi. 

Where A-Isoself c = {A;| — y} is the set of cuspidal angles, and A-Isoselfi is the set 
of the roots of P{tp)- 

4.1.2 Characterization of Isoself c on generic points of C. 

We have to solve the equation : 



Y^fifr + 7T + S)h* = YJii<p)h* + 0(h 8 ), 



i=4 i=4 



for small enough values of 5 and h. 

We know that this equation is equivalent to : 



fe(<P + fM + h (f 7 {<p + 5 + n)- f 7 {<p + n)) + h (f 7 (cp + tt) - f 7 {cp)) + 0{h 2 ) = 0, 

(7) 
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For simplicities sake we denote — by ' ; hence |7| writes : 

E c = = 6f t (<p) + + *]%{<p) - hSf'M - 2hf r (<p) + 0(5 4 ) + 0(h5 2 ) + 0(h 2 ) 

By (B.P.) /» = (since 9 3 = 0). Checking /™(p) A /£ (p); we obtain : 

/6"(0)A/ f ;'(0) = -583200|r3|7r 2 

Hence : 



E c = 



E c A/;"(0) = 
E c AfZ(0) = 



Expanding E c we obtain 



E e = S/ e (Q) + 6<pfi(Q) + &f?(0) + 0(V)+ 

f /e'(O) + fW^O) + 0(*V) + f/T(0) + 0(5*) - hSf' 7 (0) + 0(h5 V ) 

-2/1/7(0) - 2 /i y> / 7 (0) + 0( V) + C(/i 5 2 ) + 0(h 2 ) = 0. 
Hence : 

E c A tf(0) = f^ 2 /;"(0) A (0) + f ^"(0) A £(0) + f /;"(0) A (0) 

-2/i / 7 (0) A /g (0) + O(M) + 0(/i <p) + 0( V) + ( J4 ) + °& 2 $ 2 ) + °( h2 ) = 0- 
and, 

^A/;"(o) = (^ + f)/;'(o)A/;"(o)+ 

o(/i) + o(V) + °(^ 3 ) + o(<5V) = 0. 

We can solve [8] in /t (implicit function theorem). 

We claim that / 7 (0) A fl (0) ^ 0. 
In fact : 

^(/ 8 A/ r )(0) = /J(0)A/ 7 (0), 
And zero is not a root of P(y) by theorem |37 
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(9) 



Replacing h by its value in ^| we obtain : 

¥>=-! + 0(5% 



t _ 1 / 6 (0)A/ 6 (0) 5 3 , ^4N 



Remark 4.2 Since h has to be positive, only one half of the curves defined by [7J must 
be considered : either 5 > or 5 < 0. 



Now let us show that this self-intersection is transversal. For this, the following ex- 
pression (the transversality-rate) T(ip, h, 8) has to be nonzero on the self-intersection. 

T(<p, h, 5) = (/ e ( v + 5) + h f 7 (<p + 5 + tt)) A (f'M + h f' 7 (tp) + 0(h 2 )). (12) 
Taking into account |IT], and expanding [T2|, we obtain : 



Since, 



T(<P> h, S) = ~tf{0) A fe(0) + 0(5') (13) 



■^"(0)A^(0) = 2/ l A( / ; A/7 ) (0 )^0 



Hence 113" 



is nonzero. 



4.1.3 Characterization of Isoselfi on generic points of C. 

Now we assume that ^ ^ 0. 
We want to solve again : 

f 6 (<p + n + 5) + f 7 (ip + 7r + 6)h = f 6 (<p) + h f 7 (ip) + 0(h 2 ), (14) 

for tp close to zero, zero being a simple root of P((p), 5 and h small. 



Expanding [TJ], we obtain : 



= 5 fM + 0(6 2 ) - h f 7 (cp + 5)-h f 7 (<p) + 0(h 2 ) 

o = 5 fM - 2 h Mv>) - h ( Mv> + 5)- Mr)) + °( h2 ) + <>(&). 

This last expression rewrites : 

E 1 = = 5 f' 6 (0) + 5 ip f'^0) + 0( V) " 2 h f 7 (0) -2hcp f' 7 (0) + 0(h 2 )+ 
0{5 2 ) + 0{hp 2 )+0{h5), 



(15) 



(16) 



19 



Checking f' 6 (ip) A / 6 (</?) we obtain : 

f'M A fl M = 32400 |r 3 | 2 vr 2 sin 2 (3 <p + 9 3 
Since #3 has to be nonzero (far from cusp), 



/ 6 (0)A/ 6 (0)^0. 



Therefore : 



E 1 = 



£7iA/b(0) = 

£?! a ;;(o) = 



Hence : 



= Ex A = J/ 8 (0) A #(0) - 2h/ 7 (0) A tf(0) + 0(h 2 ) 
+0(5 2 ) + 0(M + O(hS) + 0(5ip) 



and 



(17) 



— Ex A f' 6 — StpfiQS) A £(0) + f (0) A /;(0) 

-2/^^(0) A /^0) - 2 V/;(0) A £(0) - h5f 7 (0) A /;(0)+ (18) 
0(5 2 ip) + 0(5 3 ) + 0( V) + 5 2 ) + 0(5 <^ 2 ) + 0(/i 2 5) + 0(/i 3 ) 

If / 7 (0) A ft(0) ^ 0, (i.e. / 7 (0) ^ and / 7 (0) = A/ 6 (0)), then gives us : 

5 = 2 /7nl A^7nl * + h 0l &> V = 2X h + h 0\<P, h), A ? 
Replacing S by its value in [l^ and dividing by h (h > 0) ; Ex A /g(0) = become : 

2 % (/ e A /r)(0) " 2 /r(0) A ^ (0) + 2A h T {fe A /7)(0) + 0(/ " 2) + 0{ip2) = (19) 

Recalling that / 7 (0) A f' 6 (0) = 0, and solving [19| in tp\ we obtain at the end : 

(p= -\h + hO(h) 

(20) 

5 = 2\h + hO(h) 

Since h > the curve defined by |2^ is a smooth one starting from zero. 

We have to show again that this self-intersection is transversal. 
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The transversality-rate T ± ((p, h, 6) has to be nonzero on the self-intersection. 
TM h, 5) = (f' 6 (tp + 5) + hf' 7 (v + 5 + tt)) A (f'M + h f'M + 0(h 2 )). 



TM h, S) = (f'M - h f'M +8fi(tp) + 0(5 2 ) + 0(6 h)) A (f'M + h f'M + 0(h 2 )). 
It is easy to see that : 

TM h,S) = h ^(/ 6 A f 7 ) + 0(/* 2 ) = 0(h 2 ). (21) 

Since zero is not a double root of f' 6 A /V, then T(<p, h, 8) ^ for ^ close to zero. 
From now we assume that /7(C)) = 0. 
Hence : 

^(/6A/r)(0) = /,(0)A/ 7 (0)^0. 

In this case we will need higher order jet of CL : namely /g(y). 
Remark 4.3 / 8 insufficient because of the fact that f 8 (tp + n) = /s (<£>)• 
We have to solve now : 

K = 5 fM + f fM + f + o(5 4 ) - 2/1/7^) - hsfM+ 

h 2 Sf' s (0) + 0(/i 3 5) + 0(hS 2 ) - 2h 3 f 9 (0) + 0(h 4 ) = 0. 

/e(0) A /r(0) 7^ 0. 

Hence : 



£?! = 

Also : 



a fM = 
a / 6 (o) = 



and 



= E[ A /; = 5 f' 6 (0) A f'M - 2 / 9 (0) A / 7 (0) + 0(/i 4 ) + 0(5 2 ) + 
0(iph) + 0(hS) + 0(Sip). 

= E[ A f' 6 = -2hcpf 7 A /; - 2 /i 3 / 9 (0) A ^(0) 
O(V) + 0(5 2 ) + 0(<5p) + 0(/i<5) + 0(h A ) 
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(22) 



(23) 



(24) 



From |23| one obtains : 



S = 2 /9(0)A/ 7 (0) h * + h * oH h) 



Replacing 6 by its value in |4], and dividing by h (h > 0), E[ A /g(0) = become : 
-2p/ 7 (0) A /i(0) - 2 /i 2 / 9 (0) A / 6 (0) + 0(h 3 ) + 0(^ 2 ) = 

Hence : 

cp = A'h 2 + h0 2 (h) 

(25) 

5 = 5' h 3 + h0 3 (h) 
Again the curve defined by ^ is a smooth one starting from zero. 

Showing now that this self-intersection is transversal. 

The transversality-rate T^ip, h, 5) has to be nonzero on the self-intersection. 
T 2 (<p, h, 5) = (f' 6 (ip + 5) + hf' 7 (v + 5 + tt)) A (f'M + h f 7 {<p) + 0(h 2 )). 



T 2 (<p, h, 5) = (f'M - h f 7 (y) + 0(5) + 0(6 h)) A (f'M +hf r (<p) + 0(h 2 )). 



As in the previous case : 

r 2 (p, h,S) = h ^-(/; A f 7 ) + 0(h 2 ) = 0(h 2 ). (26) 

Zero is not a double root of / 6 A hence T 2 (ip,h,5) ^ close to zero. 

4.1.4 Characterization de Isoselfi on isolated points of C 

Double roots of P(y?). 

We recall hypothesis : 

i) iif'e A / 7 )(0) = (0) A / 7 (0) + /£(0) A / 7 (0) = 0. 

ii) fl 3 ^ 0. 

The equation of the self-intersection is the following : 

El = 5f'M + f (<p) + § + 0(5 4 ) - 2hf r (<p) - 

(27) 

/i 2 5/s(0) + 0(h 3 5) + 0(M 2 ) - 2h 3 f 9 (0) + 0(h 4 ) = 
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Our hypothesis ensure us that : 

/ 8 (0)A/?(0)^0. 

Hence : 



E{ = 



We interesse at first to E\ A / 6 ' (0) = 0. 
We obtain : 



£w;'(o) = o 

E\ A / 6 (0) = 



= ^(0)A/;'(0)-2/ i / 7 (0)A/;'(0) + 

0{5 2 ) + 0{5ip) + 0(h 2 ) + 0(h(p) + 0(h5) 
Solving this last equation in 6 ; one obtains : 



'"m^ h+ ^ (28) 

We will distinguish two cases. 

i) Mo) ^ o. 

In that case : 

/ 7 (0) = A/ 6 (0) with A ^0. 

and (29) 
5 = 2\h + 0(h 2 ) + 0(hip). 

We will consider now El A / 6 (0) = 0. Assume that there is at least one solution. 
This solution is on the form: 

V = x h + X 2h 2 + X2 h 3 + 0(h 4 ). 

It is easy to see that x 0- 

Replacing ip by its value and dividing by h? we obtain : 

m x 2 + 2Xmx + n + O(h) = 0. 
where, 

m = A^'(0)A/ 6 (0) + / 6 (0)A/ 7 '(0) 

n = f fi{0) A / 6 (0) - 2A£(0) A / 6 (0) - 2A/ B (0) A / 6 (0) 
By remark RT3I we have the following important fact : n ^ 0. 
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Hence : m^O. 

Denoting xi an d X2 two presumed solutions of 

m x 2 + 2Xmx + n = 0. (30) 

Then : 

Xi + X2 = ~2X (31) 



If Xi — X2 = Xo — — A, then —m\ + n = 0. This is impossible by remark |37| once 
more. 



Therefore has not a double root on M. 

Hence in this case there is either no self-intersection or two self-intersections which 
satisfy the following equations : 

5 = 2Xh + 0{h 2 ) 

tp = xh + 0(h 2 ) where, x = Xi e ^ £ = 1,2 

It remains to show that theses self-intersections are transversal. 
Denoting by T^(<p, h, 5) the transversality-rate, we have as usually 

T?(<p, h, 6) = (fa + 5) + hf' 7 ( v + 5 + tt)) A {f'M + h f'M + 0(h 2 )). 

Expanding 7i(<p, h, 5), we obtain : 

7?{<p, h, 5) = 2h\ x + A)^(/ 8 A / 7 )(0) + 0(h 3 ). (32) 

Since (% + A) 7^ 0, then 1x~((p, h, 5) ^ on our two self-intersections, hence they are 
well transversal. 

2) / 7 (0) = 

Now A = and f' 6 (0) A f' 7 (0) = 0. 
Set f 7 (0)= 7/^(0) 

Remark 4.4 Since zero is not a double root of P ((f), then 7 7^ 0. 
Taking again E| A / 6 (0),i.e. : 

= 6f 6 (0) A tf(0) - 2%/;(0) A tf(0)+ 
0(5 2 ) + O(V) + ^(/i 3 ) + O(V) + O(M). 

We obtain : 
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5 = 2-yhip + (hip 2 ) + 0(h 3 ) 
Replacing once more 5 by its value in E\ A f' 6 (0) = 0,with the hypothesis : 

V = X h + 0(h 2 ). 

Dividing by h 3 we obtain : 

m X 2 - 2/ 9 (0) A /;(0) + 0(h) = 
where 

m = 2/;'(0) A f' 7 (0) + / 6 (0) A #(0) = £(/ 6 A / 7 )(0) ^ 0. 
Thus ; X = ±v^ where re = 2 *>(°W>> ^ 0. 

Therefore if k < there is no self-intersection and if k > 0, there two self-intersections 
which satisfy the following equations : 

S = 2 1X h 2 + 0(h 2 ) 
ip = ± X h + 0(h 2 ) 

As usually we will show that these two self-intersections are transversal . 
Denoting by T 2 (ip, h, S) the tranversality-rate ; we have : 

7?(<p, h, 5) = (/ 6 ( V + 5) + h f 7 (<p + 5 + tt)) A (f'M + h f' 7 (y) + 0(h 2 )). 

Expanding this last expression, we obtain : 

7i>, h, 5) = ±2xh 2 £- 2 (f 6 A / 7 )(0) + 0(h 3 ). (33) 



Taking account our hypothesis, we know that : 

d 2 
dip 



2xh 2 ^(f' 6 Af 7 )(0)^0. 



Collision between isoselfi and isoself c . 
Now, zero is a cusp point and a simple root of P(<p). 
Hence 



/ 6 (0) = et ^(f' 6 A f 7 )(0) = / 6 » A f 7 (0) = 0. 



In compensation 
d 2 



(J , 2 (/ 6 A / 7 )(0) = f % (0) A / 7 (0) + 2/ 6 (0) A / 7 (0) ^ 0. (34) 
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and 



U (o) a f & (o) * o. 



Again, we will distinguish two cases. 




O{5 4 ) + O{h5 2 )+O{h 2 ) = 0. 



E 




A/e (0)=0 



A/ 6 (0) = 



El = Sf a (0) + 8<pfi(0) + f^"(0) + 0(V)+ 

f /e (0) + f ¥>/T(0) + 0(5V 2 ) + £ /e"(0) + 0(5 4 ) - h5 f'M + O(M^) 

-2hf 7 (0) -2 hip f 7 {0) + 0( V) + 0{h 5 2 ) + 0{h 2 ) = 0. 
/e(0) = 0, Therefore : 

^ a yj '(o) = f<^"(o) a ;;'(o) + ^f';(o) a tf(o) + f /ST(o) a yf(o) 

-2/i / 7 (0) A /J(0) - 2% + |)/ 7 (0) A #(0)+ 
0(/i <p 2 ) + 0( V) + 0(5 4 ) + <% 2 <5 2 ) + 0(h 2 ) = 0, 



Taking account the previous characterization of isoself c and isoselfi, one can 



and 



A/6"(0) = (^+?-2A/ i )/;'(0)A^"(0) 



0(M + O(V) + 0(6 3 ) + 0(h 2 ) + O(M) = 0. 



We can solve |3"7|in h. 



h = + |) + O(V) + 0(5 3 ) 



that : 
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Remark 4.5 If tp ^ =f ; then h = h(5 2 ) + 0(6 3 ). It is easy to see that this case is not 
generic. 

Replacing ip by ^ in 38 we obtain h = h(8 3 ). 
The equation KB] become : 



^/;'(0)A/;"(0) + O 3 (5) = 0. 

Which is not possible because / 6 (0) A /g"(0) ^ 0. 
Hence in that case there is no self-intersection. 



(39) 



4.2 The stability 

In this case the suspension of the exponential mapping is not a "Whitney map". The 
main argument here is the Mather theorem. 

If 02,2 = 0, then j3^s (or r^) and fi are nonzero. The intersection CL^ between CL 
and the planes w = c (for small enough values of \w\ ) is a closed curve having fold-points, 
cusp-points (six) and self-intersections (see figures). 

According to our previous subsection "characterization of self-intersection on C" (Which 
finds here its justification), all the generic self-intersections of CL^ are transversal. 

Thus CL is in "general position ", which ends the proofs of theorems |4.2j, |4.3| and [O. 



Figures. 




CL 7 : when r2 = and r% ^ 
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CL^ (an example with 7-isoself) 




CL^ (an example with 5-isoself) 
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5 Appendix 

5.1 Figures: wave front and conjugate loci 



The waves front when r 2 7^ 




Zoom on the zone C 
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The waves front when r 2 = and r 3 7^ 0. 




Zoom on the zone C. 
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5.2 Computing the exponential mapping 

Geodesies are trajectories of the Hamiltonian vector field H associated to the Hamiltonian 
H(ip) on T*M (tt : T*M -> M) 

H^) = \^{FY + ^{GY). 

The metric (F, G) is in normal form (coordinates in A4 : £ = (cc, «;) = (2, w)) 

^ = ( 1 + ?/^)^-^/3i + !(7)£, 

G=(l + ^ / 3)^-xy/?£-f( 7 )^, 
/?(0,0) =6o 

7 = (1 + (X 2 + U 2 )f3) /o x+tf^Xl^pitxty) 

Coordinates in the cotangent bundle are (p, q, r) = (£, r). (x, y, w, p, q,r) have weight 
1, 1, 2, —1, — 1, —2 respectively. 

^- 1 (i)n 7 r- 1 (0) = {p 2 (0) + g 2 (0) = l}. 

~ ~ 1 p q 

We set p(0) = cos(/?, g(0) = siny?. For r(0) 7^ 0, we set p = —r—r, p — -, q — -, ( = 

r(0) r r 

(p, g). s denotes the arclength and t the new time : dt = r(s) ds. p, q have weight 1. 
One has, for r(0) 7^ and s small : 

1- = — ~ r=i, and — (-) = -r(s) — r=1 (40) 

<is dC as r cfe 

Or, 

* ~~ ac |r=1 ' dt ~~ dz |r=1 ^ ' 



For all (f41|) can be rewritten : 

d(z, C) 



A(z, + F ^ C, «0 + O fc+1 (^, C, ™), (42) 



i=3 



where is homogeneous of degree i, where A is a linear operator (corresponding to 
the Heisenberg sub-Riemannian metric), and where O k+1 (z, £, u;) has order (k + 1) with 
respect to the gradation : 

deg x = deg y = degp = deg q — 1, deg to = 2. 
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This can be rewritten : 

dw 
~dt 



G 2 (z, + G ^ £ w ) + O k+1 (z, C, w), 



(44) 



i=4 



where Gi are homogeneous of degree i, where G2 corresponds to the Heisenberg sub- 
Riemannian metric, and where O h+1 (z, £, w) has order (k + 1) w.r.t. the gradation. 
Initial conditions are 



Therefore, 



z(0) = 0, w(0) = 0, C(0) = (p cos(p, p simp). 



(z, Q = p ( Zl (t, p), (i(t, p)) + E P l P), Ut, P)) + 0(p k+1 ), 

i=3 



w = p 2 w 2 (t, p) + Y.p 1 Mt, <p) + o( P k+1 ). 

i=4 



(45) 



(46) 



( Zl , Ci)(t, </>) = e At (*(0), C(0)), w 2 (t,p) = / G 2 (^(r,v9), d(r,p))dr. (47) 

./o 

These last expressions can be easily computed. They give the exponential mapping for 
the Heisenberg metric : 



' zi(t, p) = 2 sin(f) (cos(y? - §), sin(p - §)) , 
Ci(t, ^) = cos(f) (cos(y? - §), sin(c/? - f )) , 
-U7 2 (t) = l(t - sint). 
Also F 3 and G 4 don't depend on to, hence, setting A = (z, (), 

A 3 (t,p) = (z 3 ,( 3 )(t,<p) =J*e A ( t -T)F 3 (A 1 (T,p))dr, 
Ml V>) = H {^(Mr, P)) ■ A 3 (r, p) + G 4 (A 1 (r, p))) dr. 



(48) 



(49) 



The following terms are easily computed, on the same way. We give the expressions that 
we shall need, F 4 , G§ don't depend on w, 



A 4 (t, p) = (z 4 , Q(t, p) = £ e A ^) F4 ( Al ( T) V9)) dT> 
U7 6 (f , = £ (^(Ai(r, ?)) A 4 (r, p) + G 5 (A 1 (r, p))) rfr. 



(50) 
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' A 5 (f , <p) = ft e A ^) (F 5 (A 1 (t, p), u*(t)) + ^(At(r, V>)) Mr, p)) dr, 
< w 6 (t, ^ = 1*1^. (A 3 (r, p), A 3 (r, p)) + G 6 (A 1 (r, p), w 2 (r))+ (51) 
^(A^r, p)) . A 5 (r, + ^(A^r, p)) . A 3 (r, p))dr. 

' A 6 (^) = / *e^(^)(F 6 (A 1 (r^),«, 2 (r)) + ^(A 1 (r^))A 4 (r, p)+ 

^(A 1 (r^))A 3 (T^))dT, 

< A 7 (t, p) = £ (FtCAxCt, p), toaCr)) + ^(^(r, p)) A 5 (r, p)+ (52) 

^(A 1 (r^))A 4 (r,^) + 

k k&(M(r, <p)) (A 3 (r, p), A 3 (r, p)) + ^(Ai(r, <p)) A 3 (r, ^))dr, 
As we shall see, we will need to compute these values for t = 2 tv only. 

5.3 The exponential mapping in suspended form, the conjugate 
loci 



z(t,p, p) =pz 1 (t, ip) + J2p l z { (t, p) + 0(p 8 ), 



S{t,p) : < 



i=3 



w(t, p, tp) = p 2 w 2 (t) + E P i w{t, p) + 0(p 7 ). 

i=4 



Integral expressions of 2j, to, were computed in our appendix 5.2 



Let us consider the variable e = ±1, according to w > or w < 0. Let us set 
t = 2ne + r (the conjugate new time has an expansion t c = 2tt e + 0(p 2 ), as is shown in 
our previous papers |J, and this will appear again here in). Therefore, r c , the conjugate 
time r, has order p 2 . 

w 2 (t) = ~{t - sin*). w 2 (t) = w 2 (r) = ne + ~(r - sinr), 

hence, wj(0) = vre, wj'(0) = 0, u^' (0) = 0. Therefore, we obtain the following important 
fact : 



w 2 (t) = 7T e + 0(r 3 ) = u£(0) + r 3 ^(r), 
for some smooth function ijj. 



(53) 
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Let us set : 



fc=,/^- (54) 



w{r, p, <p) = vrep 2 (l + 0(r 3 ) + ^(p 2 w 4 (r, <p) + p 3 w 5 {t, tp) + p 4 iw 6 (r, <p) + 
p 5 u, 7 (r, <p) + p 6 w 8 (r, ^) + 0(p 7 )). 



/i = p (1 + 0(r 3 ) + — (p 2 w 4 (r, <p) + p 3 w 5 (t, <p) + p 4 w 6 (r, <p) + 0(p 5 )))5 . 
7T e 

Using the implicit function theorem, we can solve this last equation in p. After straight- 
forward computations, we obtain : 



0(r 3 ) + 0(/i 5 ) + 0(/i 2 r 2 ) + 0(/rV) ). 



c?r fc ' * dr k 
We also denote z\ for etc. . . . 
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(55) 



Let us set now 

r = 9h 2 + ah 3 , (56) 

for some constant 9. 

p = h(l- hh^ - h 3 e^ + h\^f - - <'(0)f? )+ 
c 5 /i 5 + c 6 /i 6 + 0(/i 7 )). 

where : (57) 
C5 = — e -|^+ terms depending on the invariants fa, (i < 2) 
c 6 = — e + terms depending on the invariants (z < 3) 

Otherwise, 

9 

z(t, p, <p)=p Zl (T, <p) +J2 z i( T i f) + °(p 8 )- (58) 

We set 

(fc) d fc £i (fc) d k Wi 
A = -rr, w ; 



Replacing fl56|) , (|57f) , in this expression (p8|) , we obtain, after tedious computations : 



z{a,h,<p) =J2 Ati + Oih 10 ), 



i=3 



^3 = {6 z[ + Z 3 )\ t=27T e, M = (az' 1 + Z 4 )\t=2ne, 

A 5 = (^(9w 4 z' 1 + 3w 4 z 3 ) + \z[ +9z' 3 + z 5 )\ t =2ne 
A e = ^(-£(j w 4 z x — e 9ws z[ + 2na9z' 1 — 3e w§z 3 + 2 iraz 3 — 
4ew 4 z 4 + 2n9z' 4 + 2irz 6 )\ t=27Te ^aA\ + A% 
A 7 = a 2 A 2 + aA 7 + A° 7 , 

A 7 = ) U=2tt, A\ = i(eW 5 ^ + 27TZ' 4 )\ t =2n 

A° 7 =(±(£w 2 -e9 2 w' 4 -e9w 6 )z 1 - \9 2 w 4 z[) + 3 tt 9* z" + 
Y^w\z 3 + 3(-^w 4 + e9w' 4 - ew 6 ) z 3 + 3 e 9 w 4 z 3 + tt 9 2 z' z - 
Aew 5 z 4 + 5ew 4 w 5 + 2tt 9 z' 5 + 2tt z 7 )\ t=27V - 
A 8 = aAl + A° 8 , A 9 = a 2 A 2 + aAl + A° 9 . 

As we know from ||, the conjugate time r c is obtained for 

9 = — 6 tt e b 
To see this it is sufficient to compute 9 for 

= 4'(2 7 r6 ) ^)A(0M^) + M^ )) 

1, z\ (2tt6, tp) A nziizmuLi. — 6 7r e &o- 



9ip 



but, (27re, 99) A 
Also, we have : 



<9<£> ' dip 

-6 7r e 6 4' (27re, 99) + z\ (2ire, p) = 



(59) 



(60) 



(61) 



as it is easily checked. 

Therefore, we have the expression of the exponential mapping, in suspended form, in 
a certain neighborhood U of the conjugate locus at the source, for h (or p) small enough 



z{a, h,p) = J2 Ah' + Oih 10 ), 



(62) 



i=4 
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where are given in fl60l) . 

We have to compute the expression of the conjugate time a c in terms of h and <p. For 
this, we have to solve the following equation in a : 



dz dz 

— (a, h, <p) A — (a, h, ip) = 

Solving this equation in a, with the implicit function theorem gives : 



(63) 



c" 



i=0 



Remark 5.1 o~\ = — (z^(27re) A 9 -^ m )+ terms depending on the invariant (3{, (i < 4) and 

a\ = — (z^(27re) A d -§^ l ) + terms depending on the invariant (i < 5). We don't need the 
complete expression of those two coefficients. 

Now we give the expression of o l c {% = 1, 2, 3). 

" a° c = -(4(2-) A fg), a\ = -&(2*e) A 9 -0) 

a 2 = (-(z'^e) A d -§)- (Al A fj(27re))+ 

(4 A %){2ne){A\ A ^(2vr e ) + z'^e) A $)) 

a c 3 = -(z'^ne) A d -§)- (A 61 A ^) - (A 71 A f^vre))- 



(64) 



(4 A ^) 2 (2ne) (2 A? A ^ + £(2ice) A 
(4(2vre) Af)KA^ + 4(2-) A + 

(4 A 9 ^){2*e){A\ A ^ + z[(2ne) A M + 2 A* A ^1) 

It remains to replace the expression of a c in the expression of z in (52), to obtain the 
expansion of the conjugate locus : 

9 
i=4 



The expressions of the /, (for z < 4 < 6 and e = 1) are given in \2.2\ If we denote by 
ff the /j for e = — 1, we obtain : 
ff = -f (for 4 < i < 6) 

The expressions of f-j can be computed in the same way, just replacing flBlp in 
This has been done with Mathematica. Let : 

Pa = £44 {x 2 + y 2 ) 2 + a 44 (x 4 + — 6 x 2 y 2 ) + 4 6 44 xy (x 2 — y 2 ) + 
c 44 (x 4 - y 4 ) - 2 d 44 x j/ (x 2 + t/ 2 ) . 
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Assuming r 2 = 0, we obtain : 

1) h — 3 7r(— 21 c 44 cos(y?) + 35 a 44 cos(3 tp) — 7c 44 cos(3 y?) + 21 a 44 cos(5 </?) + 
3 r\ cos(p — 2 9i) + r 2 cos(3 99 — 2 #1) + 21 a? 44 sin(^) — 35 6 44 sin(3 <p) + 

7d 44 sm{3<p) - 216 44 sin(5 - 12 ixr\ sin(<£ - 2 0j) - 12 yrr 2 sin(3 y? - 2 #1)), 

3 7r(21 d 44 cos(9?) — 35 6 44 cos(3 p) — 7 a? 44 cos(3 + 21 6 44 cos(5 ip)— (65) 

12 7r r 2 cos(p — 2 6i) + \2t\t\ cos(3 95 — 2 #i) + 21 c 44 sin(y?) — 35 a 44 sin(3 p) — 

7 c 44 sin(3 p) + 21 a 44 sin(5 p) — 3rf sin(p — 2 6i) + r\ sin(3 tp — 2 ^)) 

2) fr + f7 = -144 7T 2 sin(2(v9-^))(cos((^),sin(v9)) = d 7 . 

It is not necessary to compute the expression of / 8 , the only thing we need is to know 
that 

f 8 (p + n) = f 8 (p). (66) 
As we know, this is stated in 0, and we verify it by Mathematica. 

For f 9 the term z 9 appears in Ag (formula (|60|)), but also, it comes through the term 
A\ /i 4 of z{a, h, p) in : A\ — (a z[) (2 tt e). 

Therefore,the expression of fg contains also the term 

T = a c 5 4(2vre) h 9 . 



This term will not play any role : the only thing that we need is to prove the remark |3T3 . 

-q— = —180 |r 3 | % sin(3y9 + ^3) (cosp, smcp), z x ^2ixe) = (cosp, sin^). (67) 

Hence, — — A fg(p) = -7— A z 9 (2ne) + other terms that depend on the invariants /3,, 
op op 

(i<5). 

We have : 

^1 A Z7{27Te ) = ^ A f 2 - e Mt~r) F (A ( j W2{t)) dT+ 

dp dp (6g) 

terms not depending on (3 G 
This last expression has been computed by Mathematica, and we don't give it here. 



We will give now the expression of the polynomial P(p) in the section \2~2 
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9 'fa 

Taking in account B7I and KoT we see that — — A dy = 0. Therefor 

dip 

Hence (up to sign) the polynomial P(ip) does not depend on e. Using Mathematica again 
we obtain : 



P((p) = — 7rf 44 cos(2 ip) + 7b 44 cos(4c/?) — 7c 44 sin(2 <p) + 7a 44 sm(4<y?) + 
r{ sin(2(^-^)). 
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